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Abstract 

We compute the stationary profiles of a coherent beam of Bose-Einstein condensed atoms propa- 
gating through a guide. Special emphasis is put on the effect of a disturbing obstacle present in the 
trajectory of the beam. The obstacle considered (such as a bend in the guide, or a laser field per- 
pendicular to the beam) results in a repulsive or an attractive potential acting on the condensate. 
Different behaviors are observed when the beam velocity (with respect to the speed of sound), the 
size of the obstacle (relative to the healing length) and the intensity and sign of the potential are 
varied. The existence of bound states of the condensate is also considered. 
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1 Introduction 



The field of Bose-Einstein condensation (BEC) of atomic vapors is undergoing rapid experimental 
development providing a rich new phenomenology and also allowing to test concepts that originated 
in other fields (mainly in the theory of superfluidity, in nonlinear optics and in atomic physics). 

Along these lines, the possibility of building atom lasers by guiding condensed particles through 
various geometries opens up the prospects of a rich variety of interference, transport and/or coher- 
ence phenomena. Cold atoms have already been propagated in various guides (see e.g. and 
references therein); more and more efficient coherent sources of atoms have recently been designed 
(using various output coupling schemes, see Refs. 1^-^) and continuous guided beams of condensed 
atoms will be accessible in the near future (see the preparatory study ||lO||). 

The purpose of the present work is to explicitly determine the different propagating modes of a 
beam of condensed atoms through a guide, as a function of the various external control parameters. 
We consider the simplest geometry, in which a guide connects two reservoirs, and treat the case of 
atoms with a repulsive effective interaction, such as rubidium and sodium. First, the transmission 
through a straight guide is reevaluated: small amplitude density oscillations, cnoidal waves, and 
gray solitons are possible propagating modes. The main part of the paper is devoted to the study 
of coherent transmission modes in the presence of an obstacle. We find that, as a function of the 
speed of the incoming beam and the size and depth of the perturbing potential, many different 
transmission modes exist. For example, some are soliton-like modes (with a peak or a trough in the 
density) that are pinned to the obstacle. They may also have density oscillations in the region of 
the obstacle. On the contrary, other modes are step-like shaped. In general, the modes may or may 
not have a wake. The wake, however, always precedes the obstacle (it never occurs down-stream). 

The obstacle is represented in our treatment by a potential that acts on the condensate along 
the axis of the guide. It can be repulsive or attractive. There are at least two simple ways to 
realize such an obstacle experimentally. The first is to bend the guide: a potential proportional 
to the square of the local curvature is created. This potential is always attractive. In view of 
future applications of the atom laser to more complicated geometries, the presence of bends seems 
unavoidable and their study is therefore far from academic. The second possibility is to illuminate 
the beam with a detuned laser field. Contrary to the case of a bend, attractive and repulsive 
potentials may be realized by varying the frequency of the laser. This latter method also has the 
advantage of allowing a better control of the relative speed between the obstacle and the beam (if 
using an acousto-optic deflector). 

In the case of attractive potentials, aside from the transmission modes there may also exist 
bound states in the guide, in which condensed atoms are trapped without possibility of escape. 
These bound states are also analyzed here. In the case of a bend we show that, due to the intra- 
atomic interactions, the bound state can support only a limited number of condensed atoms, which 
is typically of the order of 100 for rubidium and sodium. However this number can be made much 
larger for an attractive potential originating from a red detuned laser field. 

The paper is organized as follows. In Sec. 2 we set up the theoretical framework and notations. 
Our approach for describing the condensate motion is based on a one-dimensional reduction of the 
Gross-Pitaevskii equation The solutions described here can in some instances be explicitly 

written in an analytic form, for example, in terms of elliptic functions. Although we use this 
opportunity in some cases, we have chosen to keep the discussion at a qualitative level. This allows 
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us to cover a large range of experimental situations and gives a global view of all the possible 
solutions in the different regimes of the control parameters. In Section ^ we study the existence of 
a bound state of the condensate created by an attractive potential. We determine the maximum 
number of atoms the bound state can accommodate. In Sees. we study the propagation of a 
condensate through a guide connecting two reservoirs. The analysis relies on an interpretation of the 
Gross-Pitaevskii equation in terms of a fictitious classical dynamics. We focus on different specific 
examples, starting with a straight wave guide without any potential (Sec. ^). We then study the 
motion in the presence of an obstacle represented by an attractive (Sec. ^) and a repulsive (Sec. ^) 
square well. Sec. presents two alternative treatments for determining the transmission modes 
in the presence of an obstacle. The first one is a perturbative approach; the second approximates 
the external potential by a 5 function. Sec. ^ contains an analysis of the results in view of future 
experimental realizations. We present our conclusions in Sec. ^. Some technical aspects concerning 
the adiabatic approximation used for deriving the ID reduction of the Gross-Pitaevskii equation 
are included in the Appendix. 



2 An effective one-dimensional equation 

We consider BEG atoms at zero temperature confined to a wave-guide. Let x be the coordinate 
along the axis of the guide (which is possibly bent) and rj_ a perpendicular vector giving the 
transverse coordinates. We work in the adiabatic regime where the local curvature k{x) of the 
longitudinal motion of the atoms is small. If the transverse extension of the wavefunction is denoted 
by R±, this is more precisely defined by the limit R±k <C 1 and \dK/dx\ <^ 1 (see Appendix A| ). 



In this regime one can consider x and r± as Gartesian coordinates, i.e. for instance the volume 
element d^r is approximately dxd'^r±. 

It is consistent with the adiabatic approximation to make an ansatz for the condensate wave- 
function ^(r, t) of the form |11| 

^{r,t)=^{x,t)c^{f±;n) (1) 

where cj) is the equilibrium wave function for the transverse motion, normalized to unity, J d'^r±\(j)\'^ = 
1. tp(x,t) describes the longitudinal motion, and the density per unit of longitudinal length is 

liTr|2 



n{x,t) = J (iV^I^E'r = l^(a;,i)r ■ 

Notice that the transverse wave function depends parametrically on n{x,t). The adiabatic approx- 
imation is in fact a local density approximation in the sense that one assumes that the transverse 
motion is not affected by densities at points other than x. As noted in Ref. [0] this corresponds 
to the assumption that the transverse scale of variation of the profile is much smaller than the 
longitudinal one (and corresponds indeed to the limit R±k <^ 1). 

The beam is confined by a transverse potential V±{f±). Keeping in mind experimental realiza- 
tions, we often consider below the particular case of a harmonic trapping V±{r±) = r'j_ {uj± is 



the pulsation of the harmonic oscillator; we set units such that h = m = 1). As shown in [Appendix 



within the adiabatic approximation the presence of a bend results in an attractive longitudinal 
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potential V|| (x) given by 

V||(x) = -k2(x)/8. (2) 

The Gross-Pitaevskii equations for the condensate are derived through a variational principle. 
One extremizes the action 



^4 



j d^rdt{^*dt'^ - ^Ot^-*) - J dtSm , (3) 



where is the energy functional 

= j d\ Ql W|2 + 27ra,e|^'|' + + ■ (4) 

ttsc is the s-wave scattering length of the inter-atomic potential (which is represented hy a, 5 func- 
tion interaction). In all the present work, we consider the case of repulsive inter-atomic interac- 
tions, asc > 0. The extremization of 5 with Lagrange multiplier e(n) imposing the normalization 
J d^r^l^'p = n for each x (more precisely for each n) leads to the following equations 

-^Vi0 + {V^ + Airascn = e(n) , (5) 



and 

-=:^■',,i^ + iV^^+e{n))^P = i^ti;. (6) 



1.2 



In d) the non-linear term e(n) (remember that n = |V'P) is determined as a function of n from 
Eq. (|5|). In the low density limit ascn *C 1 the non-linearity in (|5|) is small. In this case, a 
perturbative solution of (|5|) leads to 

e(n) = eo + 2 ascn/a]_ , (7) 

where eo is the eigen-energy of the ground-state (j)o of the transverse unperturbed Hamiltonian 
-iVi + V±, and aj^ = 2 7r f \cl)o\'^d'^r±. For a harmonic confining potential a± = oo^^"^ is known 
as the oscillator length. 

In the opposite large-density limit nagc ^ 1 the Thomas- Fermi approximation holds, namely 
the kinetic term in (|5|) can be neglected, e is obtained as a function of n through the relation 
^tf{^) = '^CLsc'n, where Ntf is the integrated Thomas-Fermi density of states A^rF(e) = ~ 
V±)@{e — V±) d'^r±/{2Tr). For a harmonic confining potential this reads 

e(n) = 2uj±^/nasc + eo • (8) 

We remark here that the Gross-Pitaevskii equation is valid in the dilute gas approximation, 
when the 3D density n^o of the gas satisfies a^^ ^ 1 [12|. This reads here nusc ^ {a±/asc)'^^'^ 
{u = 1 in the dilute regime and u = 1/2 for high densities). a±/asc being typically of order 10^, this 
condition will be considered as always fulfilled, even at high longitudinal densities, when na^c ^ 1- 



4 



On the other hand, the weakly interacting ID Bose gas picture also breaks down at very low 
densities, in the Tonks gas regime (recent references relevant to this discussion are [1^1^). This 
occurs in the regime na^c <C (asc/a±)^ ~ 10~^ which we thus discard from the present study. 

Eq. (^), together with Eqs. (Q) and (P), are the main results of this section. They provide an 
effective one-dimensional equation for the description of the dynamics of the condensate along the 
guide. Notice that in general the non-linear term in Eq. (P) does not have the standard cubic form 
of the one-dimensional Gross-Pitaevskii equation. This happens only in the low density regime 
ttscn <C 1 when the non-linear potential is proportional to 

We emphasize that Eq. (^) relies on an adiabatic approximation. It is well known that in many 
instances this approximation gives accurate results well beyond its strict domain of validity. Two 
examples relevant in the present context are Refs. and [|18|, where the propagation of waves 



(without non-linear effects) was studied in the extreme non adiabatic case of wave guides with 
a discontinuous curvature and a sudden constriction, respectively. The adiabatic approximation 
was nevertheless shown to be applicable in these systems. On the basis of these examples (and of 
others) , one can consider that the results presented here have a wide range of validity. 

For practical purposes it will appear useful in the following to introduce a longitudinal healing 
length ^ defined for a constant longitudinal density n as 

-L = e(n) - eo . (9) 

This gives ^ = ^a±{nasc)^^^'^ in the low density regime and = ^a±{nasc)^^/'^ for high densities 
in a harmonic confining potential. 

We now study in detail the different solutions of Eq. (^). Although the attractive potential 
V|| (x) appearing in the equation of motion was due to the presence of a bend in the guide, our 
results are very general and V|| (x) could be of a completely different physical origin. In particular, 
we will consider in Sec. ^ the case of a repulsive potential which cannot be produced by a bend. 

3 Bound states 

We first study the existence of bound states in the guide due to an attractive potential V|j(2;). 
The existence of bound states in the quantum mechanical motion of non-interacting particles in a 



bent wave guide has been extensively considered in the past (see e.g. [19-21| and Refs. therein). 
It has been shown by Goldstone and Jaffe that at least one bound state exists in two and three 
dimensional bent tubes (cf. the discussion in [Appendix A ). The particle is trapped because its 



energy is lower than the first propagating mode of the straight guide. In the case of a condensed 
beam, we are interested in whether bound states occur in the presence of interactions. 

In the extreme dilute limit agcn Eq. (|6|) reduces to an ordinary one dimensional Schrodinger 
equation and the existence of a bound state in V|| (x) is guaranteed by general theorems of quantum 
mechanics |2^. Hence there exists a state of the condensate whose energy is lower than the energy 
eo of the first propagating mode of the guide. This state is localized in the region where V|| is 
noticeable, and we assume for simplicity that this happens in some finite region around x = 0. 
With increasing number of atoms in the condensate the non-linear effects come into play. The 
repulsive intra-atomic interaction increases the energy of the bound state (as well as its spatial 



5 



extension), and there is a threshold beyond which this state disappears. Therefore for a sufficiently 
large number of atoms no bound-state is expected to occur. We now determine the threshold 
quantitatively by determining the maximum number of Bose condensed atoms the bound-state can 
accommodate. 

Near the threshold the state is very weakly bound, and the wave function extends over distances 
much greater than the range of the potential V||. Hence, in this limit it is legitimate to make 
the approximation V||(3;) ~ X6{x), with A = dx < 0. This approximation is not 

contradictory with the assumption of adiabaticity of the motion. 

We look at stationary solutions of Eq. @ of the form tp{x,t) = A{x) exp{—i fit), fi being the 
chemical potential and A a real function. In the regions where the potential is negligible (i.e., for 
X with the above replacement of V\\ by a 5 function) Eq. (P) can be integrated once, giving an 
equation for n{x) = A? , 



n'2 



/■n 

-|-e(n)=/xn, where e(n) = / e{p)dp. (10) 



With the convention defined in Section 2, the normalization is n{x) dx = N, where N is the 
total number of particles in the bound state. The density n can be shown to be an even function 
of X, and the matching condition at x = reads n'(0+) = 2An(0) (' denotes d/dx). Using (|lO|) 
and these two conditions we arrive at a set of two equations determining n(0) and p as functions 
of A and N, 

yn(0)=e(n(0))-Mn(0), (11) 



and 

(n/2)V2dn 



y^e{n) - fi 



n 



N . (12) 



For solving this system, one needs to know the explicit form of the function e(n). We will see that 
for realistic values of A corresponding to an attraction issued from a bend in the guide the bound 
state can accommodate only a small number of particles. In the case of a bend, it is thus sensible 
to concentrate on the low density limit (see the estimate at the end of this section). Using (^), in 
the low density regime we have e{n) = eon + ascn'^/o.^ and we obtain 

^ = ^o-2[j V^\{x)dx + ^Nj . (13) 

The first, negative, term inside the large parentheses is due to the attractive potential produced by 
the bend (cf Eq. @), while the second one comes from intra-atomic repulsive interactions in the 
condensate. Eq. (|13|) clearly displays the existence of a threshold. When the number of atoms 
occupying the bound state increases, the chemical potential fi increases and eventually reaches the 
threshold eo at which the state disappears. This occurs for a number of atoms N^ax given by 



N -< 



+00 

V|j(x) dx 

— oo 



(14) 
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For an arbitrary potential, as — > Nmax the spatial extension of the bound-state diverges. 
Thus the approximation of Vn by a (5-function' is well justified in that limit, and we expect Eq. (14) 
to be very accurate . 

The order of magnitude of N^ax can be estimated by considering a bend of constant radius of 
curvature and bending angle 0. From (|2|) and ( [l^ ) we get N^ax = o?\Q / {SascRc)- For a guide 
with i?c = 5a_L, = 7r/2 and a±_ ranging from 1 /im to 10 /im, N^ax ranges from 7 to 70 atoms 
for a condensate of ^''Rb atoms (a^c = 5.77 nm). For ^^Na, Nmax is doubled (since the s-wave 
scattering length of ^^Na is age = 2.75 nm). 

If the attractive potential originates from a red detuned laser beam, using the estimate of Sec. ^ 
(Eq. ([5^)) one obtains a± ^j(a;) dx 10® (to be compared to the value 10~^ that applies to 
a bend). From this and from Eq. (^) it follows that for rubidium Nj^^^x can be as large as 10^. 
In this case, however, one does not remain in the low density regime where Eq. (|4|) holds (this 
regime is valid if n(0) a<jc <C 1, which from (|ll]) and ( p^ gives Nmax ^ a±/a 

SC ^ 

10^). Working m 

the high density regime instead, one obtains 

3 ai 



max — „ 

lb a. 



+ 00 

V|| (x) dx 



(15) 



One gets from this equation a value of N^ax of the order of 10^*^. Note, however, that in this 
regime the high density approximation is valid at x = 0, but violated for large x (when the density 
tends to zero). Hence, without giving a precise order of magnitude, it is nevertheless clear from 
the previous estimates that the maximum number of atoms the bound state can accommodate is 
very large in this case. 



4 Transmission modes 



From now on we concentrate on the stationary states of a beam of condensed atoms connecting 
two reservoirs. For that purpose we look at the stationary solutions (in the reference frame of the 
laboratory) of Eq. (0) with ^ having a finite value at x — > ±00. We write 



tpixjt) = exp{—i fit} A{x) ex.p{iip{x)} , 



(16) 



with A and (p real functions. Since the wave function extends to infinity, the chemical potential 
satisfies fi > cq. The density is n = and the beam velocity is v = cp' . After factorizing out the 
phases, Eq. (|6|) splits into two real equations corresponding to its imaginary and real parts. The 
former imposes flux conservation, namely the product n{x) v{x) is a constant that we denote by 

$ = n{x) v{x) . 

The real part gives a Schrodinger-like equation for A{x) 



1 



A" + 



e(n) + V|| + 



$2 
2^ 



A = fiA. 



(17) 



In this section we consider the transmission modes of a condensate through a straight wave 
guide with no obstacle. This corresponds to solving Eq. ( |l7|) with Vn = 0. The modifications in the 
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beam density and phase produced by the presence of an obstacle are considered in the following 
sections. When V|| = 0, Eq. (|T7|) can be integrated once, yielding 

1 o $2 

-A'^ + W{n) = Eel with W{n) = -e{n) + fin + — . (18) 

e(n) in (^) is the same as in Eq. ( p!o[ ) and Ed is an integration constant. This constant is denoted as 
a "classical energy" because (|^) has an interpretation in terms of classical dynamics. It expresses 
the energy conservation of a one-dimensional Hamiltonian system for a fictitious classical particle 
with "position" A and "time" x moving in a potential W{n = A^), Ed being the total energy of 
the particle. The solutions A{x) therefore coincide with the "classical" solutions in the potential 
W{n = A^). The chemical potential fi and the flux $ flx the shape of the potential W{n), while 
Ed selects a "trajectory" A{x) in this potential. 

To clarify the physical meaning of Ed consider the linear (i.e., non-interacting) case e{n) = 0. 
Then the natural way to write the solution (|l6|) is a superposition of plane waves, 

jp{x,t) = ex.p{—i fit} [a exp{ikx) + P eyip{—ikx + i6)] , (19) 

where k is the wavevector and 9 an arbitrary phase. In terms of the two real parameters a and 
[3, the flux and the classical energy are written $ = /c (a^ — /5^) and Ed = k'^ (a^ + /?^). Ed is 
therefore a measure of the total intensity of the left and right incoming beams. Ed can be varied 
while keeping $ and fi constant (and therefore W{n) constant) by changing the amplitudes a and 
P simultaneously while preserving the difference (a^ — /3^)- 

For studying the shape of W{n) in Eq. ( p!8| ) in the presence of interactions it is customary to 
plot fi — dW/dn = e(n) + $^/(2n^) as a function of n, as represented in the top part of Fig. ||. At 
low densities (n 0) the term $^/(2n^) dominates. At high densities the intra-atomic interaction 
contained in e(n) takes over, and leads to a monotonic growth for large values of n (due to the 
repulsive interactions, e(n) is an increasing function of n). At intermediate values there is a mini- 
mum at a density denoted no. The relevant case (leading to finite densities at infinity) corresponds 
io fi > e(no) + <I>^/(2no), and is shown in Fig. |l|. Then the derivative of W{n) is zero for two 
densities ni and n2, with ni < no < n2. The corresponding plot of the potential W{n) is shown in 
the lower part of the figure. In order to have a finite density at infinity one should also impose a 
bounded motion of the fictitious classical particle, and this results in the two additional conditions 
(i) l^(ni) < Ed < W{n2)-, and (ii) n(x) < n2 for any value of x. 

The different types of solution A{x) are therefore described by the different motions a classical 
particle undergoes in a potential W{n) at the allowed energies Ed- The two simplest solutions 
correspond to the fixed points of the potential, n{x) = ni or n2, where the "classical" particle 
remains at rest. They correspond to constant density solutions. Since the densities are different 
(n2 > ni) and the fiux has the same value the velocities Vj {j = 1 or 2) of the condensed beam 
are also different, with V2 < vi. vi (resp. ^2) corresponds to a beam velocity above (resp. below) 
the speed of sound. To see this, we first note (see below) that for a condensate at rest (i.e., v = 0) 
with uniform density n{x) = n, the sound velocity c is defined by 

cHn)=n^. (20) 
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In the case of a moving condensate with uniform density n, one has a well defined velocity v = ^/n 
given by /X = e{n) + ^v"^ (see Eq. (|T7|) ). From Eqs. (|l8| ) and (|20|) , we have moreover d^iy/dn^ = 
(u^— c^(n))/n. Since d^l^/dn^ at ni (resp. is positive (resp. negative), it follows that vi > c(ni) 
(resp. V2 < c(n2)). 




min max 



Figure 1: Schematic behavior of the functions e(n)+(f>^/2 -n? (top part) and W{n) (bottom part) as a function 
of 71. rii {i ^ 1 or 2) is defined by t{ni) + /2nf = /i and is a zero of dW/dn; uq is a zero of the second 
derivative d^W/dn^. For ^ and $ given, a beam of uniform density has either a density ni (and a velocity 
higher than the sound velocity c(ni)) or a density n2 (and a velocity lower than 0(71.2)). At a given i?c/, 77,mi„ 
and Umax are the minimum and maximum values of the density oscillations (W{nmin) — W{nmax) = Ed). 

We now consider the density profile of the transmission modes in the vicinity of the constant 
solution n{x) = rii. For energies Ed slightly higher than W{ni) the stationary solutions are 
sinusoidal waves of the form 

/2 

A(x) =Ai + ^[Ed - W(ni)]^/^ cos (kx + e) , (21) 
k 

where k"^ = cPW/dA'^\^ = 4($^/?7,f — c^(ni)) and rii = A\ is the constant density. This sinusoidal 
wave is reminiscent of a sound wave, but sound waves are progressive whereas Eq. ( ^l| ) describes a 
standing wave. The structure of the sound wave (and the whole spectrum of elementary excitations) 
is better displayed by a slight modification of the procedure used so far. Instead of the stationary 
ansatz ( |T6|) one looks for solutions of the form 'il^{x, t) = exp{— i/it}A(x — nt) exp{ii?i>(x — nt)}, where 
w is a constant parameter, physically interpreted as the velocity of an arbitrary moving frame. Mass 
conservation now reads $ = n{v — u) (^> is the flux in the moving frame) and Eqs. ( 
the same form, with /i replaced by /i + ^i^/2, all functions now depending on X = x — ut and 
not merely on x. The constant solution n{x) = ni can now be given zero velocity {vi = 0) if one 
chooses u such that u + ^>/ni = (this is of importance because one wishes to study elementary 



9 



excitations in a system at rest). In this case, a perturbative treatment of ( p!8| ) for Ed near W{ni) 
again gives a solution of the form ( ^l]) with x replaced by X. This is a progressive wave depending 
on k X = k X — LOkt with uj^ = ku, and thus satisfying the Bogoliubov dispersion relation 

The long wave-length limit of ( p^ ) corresponds indeed to sound waves with a sound velocity c(ni) as 
given by Eq. (|20|). It is also possible to obtain in this way the dispersion relation of the elementary 
excitations of a beam moving at constant velocity vi (which is simply Eq. ( |2^ ) Doppler shifted). 

Note that our approach is unable to reproduce the decrease of slope of the spectrum of ele- 
mentary excitations that occurs in the high density regime, for wave vectors k of the order of the 
transverse extension Rj_ of the condensate. This effect, predicted in Refs. [^,|2^ and observed 
numerically in goes beyond the quasi-lD approach: it occurs when the excitation has a wave- 
length allowing exploration of side regions of the condensate that have lower local sound velocity. 
Hence, it cannot be reproduced by using the adiabatic ansatz (|l]). 

The structure of the stationary solutions for energies E^i close to (but lower than) W{n2) is 
totally different from the sinusoidal waves we just discussed (which exist for Ed >,W{ni)). The 
uniform solution n{x) = n2 coexists with a solitary wave corresponding in the classical analogy to 
a motion along the separatrix located at Ed = W{n2). This solitary wave has constant density 
n2 at X — > ±00, and a trough whose minimum density satisfies the condition W^rimin) = W{n2) 
(cf. Fig. |2|, top part). As the energy is lowered from Ed = W{n2), density oscillations appear 
whose amplitude decreases as Ed diminishes. These solutions are cnoidal waves (see, e.g., p7| ) 
with periodic oscillation between two values rimin and rimax, as defined in Fig. ||. As the energy Ed 
is further reduced the density profile tends continuously to the sinusoidal waves discussed above. 
This transition is illustrated in Fig. |2|, which is drawn in the high density regime for a transverse 
harmonic confining potential (i.e., with e(n) given by Eq. (P)). The same qualitative behavior is 
valid for any density regime. 

Figure |2| summarizes the possible density profiles of the transmission modes of the condensate 
along a straight guide. In the remaining sections we consider the modifications induced by the 
presence of an obstacle in the flow of the condensate. Finding the transmission modes now reduces 
to a scattering problem in which two of the "free" modes discussed in this section are matched 
by the potential representing the obstacle. The correct boundary conditions to be imposed are 
determined by the relative value of the phase velocity 

vp{k) = ^ (23) 

with respect to the group velocity 

Both functions start from the value c(n) at = and then increase monotonically, with Vg{k) > 
Vp{k) for any k > 0. For the stationary motion of an obstacle in a condensate at rest, Vp coincides 
with the velocity of the obstacle with respect to the beam, this is the condition of stationarity. The 



+ 



(22) 
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energy transferee! to the fluid propagates with a velocity Vg greater than the velocity of the obstacle 
with respect to the beam (since Vg > Vp). As a consequence, radiation conditions require that the 
wake is always located ahead of the obstacle (i.e., up-stream in a frame where the disturbance is at 



rest), with no long-range perturbation of the fluid on the down-stream side |28|. 




Figure 2: Beam density along a straight guide (high density regime) with V2/C2 = 0.7 (lengths are given 
in units of ^("-2), see Eq. (^). Top: solitary wave {Ed — W{n2)). Middle: cnoidal waves existing for 
W{ni) < Eel < W{n2). Bottom: for Ed ■^W{ni) the cnoidal wave deforms continuously to a sinusoidal 
small-amplitude wave. 



5 The attractive square well 

Having recalled the different solutions existing in a straight wave guide, we now discuss the influence 
on the transmission modes of a longitudinal potential representing a motionless obstacle placed in 
the trajectory of the beam. Specifically, we consider a potential V||(2;) that vanishes everywhere 
except in a finite region < x < a where it takes the constant value —Vq (Vq > 0). If its origin 
is the presence of a bend, the square well potential corresponds to a wave guide with a constant 



curvature over a finite length < x < a and straight elsewhere |29|. Apart from considerations 
related to its physical origin, this model potential is of interest because it allows one to understand 
in a simple case the different stationary regimes occurring also in more complex potentials. 

An important point for the determination of the transmission modes of the condensate along 
the guide is the boundary conditions. As discussed at the end of the previous section, amongst 
all the possible stationary solutions that exist in the presence of a scattering potential, the only 
physical ones are those that tend to a fiat density down-stream. Hence we consider density profiles 
tending to a fiat density at x — > —00, with n{x —00) — > noo, and with a negative velocity 
Voo = ^/rioo- The sound velocity at infinity Cqo = c(noo) will also be chosen negative in all the 
following. This corresponds to a beam incoming from the right, unperturbed far down-stream by 
the presence of V\\ , and characterized by the two parameters Vqo and n^o (or equivalently by ^ and 
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$, since fi 
units of ^ = 



= e("-oo) 
[2(e(no 



+ and $ 



60)1-1/2 (cfEq. (I)). 



j). Moreover, we will systematically express lengths in 



When V|| is a square well, Eq. (17) takes a particularly simple form. Everything happens as for 
a straight wave guide, except that iJ. in (^) is shifted to /i + Vq in the region < x < o". Hence, as 
in the case of a straight guide, one has an integral of motion, but it takes a different value in each 
portion of space. 



x<0 , 



0< X <a 



a < X . 



(25) 



W{n) in Eq. (25) is defined as in Eq. ([iSl). , E^^ and E^^ are the values of the integration constant 



in each region. Since the solution is flat far down-stream (when x — > — oo) one has E, = W{nc 



The matching at x 
leads to 



and X = a imposes continuity of the density and of its derivative. This 



E^i+VoniO) 



E^i 



E+ + Von{a). 



(26) 



Different types of solution satisfying Eqs. (25) exist, depending whether the far down-stream beam 
velocity Voo is greater or smaller than the speed of sound c(noo) = Coo- We consider these two 
different regimes separately. 



5.1 low beam velocity: foo/coo < 1 

The first type of solution we consider is rather intuitive if employing a perturbative treatment (see 
Sec. 1^. It corresponds to solutions with an increased density in the region of the potential. In the 
following, we refer to these solutions as "B-solutions" (where B stands for bump). 

The B-solutions are found by looking for solutions with a density increasing when x moves 
from x = — oo towards the origin. Since at a; — > — oo the density has the constant value rioo, and 
since we are imposing foo/coo < 1, ^oo coincides with the uniform density denoted n2 in Sec. 
and E~i = W{noo). But contrary to the solitary wave of Sec. Q in which the density decreases, 
in a type-B solution we move to the right of tIoo = n2 along the separatrix (see Fig. |3|) and the 
B-solution has a density peak instead of a trough. 

Since the boundary conditions are fixed down-stream, here and in all the following, we find it 
more convenient to integrate Eqs. ( p5|) starting far in the rear of the obstacle, i.e. from left infinity 
(remember that the beam is incident from the right). Starting from a value Uoo at x = — oo, the 
density has reached a value n(0) at x = 0. From this point on the equivalent "classical particle" 
moves in the potential 

Wo{n) = W{n) + Von (27) 

at an energy E'^^ = E~^ + Von(O). At x = o" it evolves again in W{n). Since to reach a finite density 
when X — > +oo one needs to have < W{noo) = E^^, from (|2^ ) this imposes n{a) > n(0). But, if 
the inequality is strict, one has E"^ < VF(noo) and n{a) > n(0) > Uoo- this case should be excluded 
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Figure 3: Type-B stationary solution in an attractive square well. The left plot displays the density profile. 
The two vertical solid lines indicate the location of the square well. This plot corresponds to = 2 ^ and 
Vb = 0.5 in a low-density beam at velocity Woo/coo = 0.7 (the healing length ^ is defined as in and is 
computed at density Uoo)- The right plot illustrates the behavior of the solution in the diagram [n, W{n)). 
The solid curve represents W{n) and the dashed one Wa{n). In the rear of the obstacle, the density evolves 
from TT-oo to n(Q) in the potential W{n) with a classical energy = W{noo), then from n(0) to rimaxfi and 
back to n(0) in the potential Wo{n) (with energy and finally the up-stream density goes from n(0) to 
Uoo in W{n). 



since from Fig. |^ one sees that this leads to a diverging density at right infinity (the classical particle 
escapes to infinity). Hence one should have n(0) = n{a), = i?^ and the B-solutions are even. 

These solutions exist for any type of attractive square well. For a given beam (characterized 
by Hoc and Voo) and given values of a and Vq, the value of A{0) is determined by demanding that 
the amplitude varies from ^4(0) to its maximum value A^axfi and back over a distance a. A^axfi 
is determined as a function of A{0) from the equation Wo(^) = -E^^ whose two smallest positive 
solutions are denoted Amin,o and Amax,o (in all the following we denote with an index "0" the 
quantities concerning Wq and defined as for W in Fig. |l|). We have 



a = 2 I — = ^2 



A(0) 



A' 



dA 



A(0) 



El - W,{A) 



(28) 



For sufficiently small a the only existing B-solution is the one described above. However, new 
B-solutions appear as the width a increases, because the "classical particle" before evolving back 
in W{n) has enough "time" to make one (or several) oscillations in WQ{n). The general density 
profile of the B-type increases from n^o at x = — oo up to x = 0, and has N maxima and N — 1 
minima between x = and x = a, with = 1,2,3... We denote this a Bjv-solution. Figure y 
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corresponds to a Bi solution. The behavior of a B2 solution is illustrated in Fig. ^ For an arbitrary 
N, Eq. dH) takes the form 



where Aminfi and Amax,o are, as indicated before, the two smallest positive solutions of Wo(^) = 
E^i- The width of the potential below which a given solution disappears occurs when ^(0) = 
■y/noo = Amaxfl- In this case the solution is perfectly flat for x < and x > a and is a portion 
of cnoidal wave (with oscillations) in the region < x < o". For a given Vq, this forces a to be 
larger than the value o"jv(Vo), 

/Amaxfi 
I • (30) 

■/^ -lE^- Wq{A) 




Figure 4: A B2 solution in an attractive square well. Left and right parts as in Fig. g. The plot corresponds 
to (T = 4^ and Vq^^ = 0.5 in a low density beam at velocity Woo/coo = 0.7. The down-stream density 
evolves from rioo to n(0) in the potential W{n) with a classical energy — VF(noo), then it makes a 
complete oscillation in Wo(n) starting from n(0) towards rimaxfl and ending again at n{a) — n(0). Finally 
the up-stream density goes from n{<T) to rioo in the potential W{n). 

Other types of solution, different from the B-family, exist for a beam velocity lower than the 
speed of sound. They correspond to density profiles that decrease from the down-stream asymptotic 
value rioo as x moves from —00 towards the origin, with = W{noo). Hence this type of solution 
is a portion of a solitary wave in the rear of the obstacle. We refer to these solutions as the 
D-solutions (where D stands for depressed). 
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For a D-solution, the density has a value n(0) < tIoo at x = and from there on, the equivalent 
classical particle evolves in the potential Wq. In the simplest case the particle bounces once on 
the repulsive core at the origin, namely the density further decreases until it reaches a value Ti^iYi^o 
(satisfying E^^ = Wo(nmm,o)) and then increases until x = a. Then the classical particle evolves 
in W{n) again, with and energy E'^ that has to be lower than or equal to W{noo), and this 
imposes n(o") > n(0) (cf Eq. (^)). Note that here, contrary to the case of the B-solution, the 
strict inequality is possible; it corresponds to < W{noo) and n{a) < n^, i.e., the up-stream 
solution is a cnoidal wave. In the particular case that n{a) = n(0) the solution is even. A generic 
density profile with a cnoidal wave is represented in Fig. ^(a), and the even solution is represented 
in Fig. |(b). 




Figure 5: A Di solution in an attractive square well. Left and right parts as in Fig. || Plot (a) is the 
generic case (with n{a) > n{0)) and (b) is the symmetric case where n{(j) = n{0). These plots correspond 
to cr = 2^ and Vo£,^ — 0.1 in a low density beam at velocity iioo/coo = 0.7. The right plot illustrates the 
behavior of the typical solution (such as displayed in part (a)) in the diagram {n, W{n)). The down-stream 
density evolves from rioo to n(0) in the potential W{n) with a classical energy E~i = W(noo). Then in Wq 
(with an energy E^^) from n(0) to the minimum density n^infi and back to a value n(cr) (larger than n(0)). 
Finally the up-stream density oscillates in the classical potential W{n) as a cnoidal wave. 

For a given well depth Vq the simple D-solutions of Fig. ^ (which we denote as the Di solution, 
with one minimum in the region of the potential) do not exist for all values of a. When the well 
becomes very large, a may exceed the period of the oscillation of density in the well (i.e., the "time" 
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period of the "classical particle" evolving in Wq). In that case the Di-solution disappears. The 
limiting case corresponds to n(0) = = n{a) i.e., to a flat density outside the region of the well 
having one oscillation in the region of the well. This upper limit is exactly the lower limit o"2(Vo) 
below which the B2 solution does not exist. As a consequence, when a increases from a small value 
the Di-solution disappears at o" = (72 (Vq) and becomes the B2-solution (of the type illustrated in 
Fig. 1^, which is then allowed for any larger value of a. The point is that for a B-solution, when 
changing from potential W to Wq in the (n, W) diagram, one can jump arbitrarily close to the 
separatrix of Wq, thus making the period in the region of the potential as large as desired. Hence, 
once a B-solution exists for a given a, it exists also for any larger value. For a D-solution however 
the period in Wq is limited: it takes its largest value if one enters and leaves the region of the 
potential Wq with a flat density (i.e., n(0) = Uoo = ^{a)). 

The existence of the Di solution for a < a"2(Vo) depends on the value of Vq and of the relative 
positions of the curves W{n) and Wo(n) in the {n,W) diagram. One regime is set by small values 
of Vq such that the condition rimin < ^^1,0 is satisfied, where n^m is the smallest positive solution 
of W{n) = W{noo) (see Fig. |l| in the case Ed = W{n2 = Uoa)) and ni_o is the first zero of dWo/dn 
(remember that we denote with an index "0" the quantities concerning Wq and defined as in Fig. || 
for W). In this case one can easily check that any value of a smaller than a < (T2(Vo) corresponds 
to an acceptable Bi solution. On the contrary, for larger values of Vq for which rimin > ''^1,0 is 
satisfied there is a minimum width c7*(Vo) below which Di solutions do not exist: 

al{Vo) = V2 , (31) 



where E^i = W{nmin) + ^o"-mm, ^min,o and A^axfi being solutions of Wo{A) = E^^ (one has 
-^max,o = Afnin)- The limiting case a = crl(Vo) is illustrated in Fig. ^. 
The situation for the solutions of type Di is summarized in Figure |^ 

D-solutions can oscillate in the region of the well, as B-solutions do. We denote by Djv-solution 
a type-D solution with N minima. Contrary to the case of B^v-solutions, there exists a maximum 
width for D^r-solutions to occur. When the D^f-solution disappears, it becomes a B^v+i-solution 
(exactly as discussed above for = 1). Djv solutions are not necessarily even as Bjv solutions are; 
ahead of the obstacle (for x > a) they typically consist of a cnoidal wave (see the discussion for 
N = 1 and Figure I). In this case we do not consider that the minima of the cnoidal wave outside 
the well increase the index "A" in the name "Djv-solution" (for instance the profile of Fig. ^(a) 
corresponds to a Di solution although an infinity of minima occur up-stream). Furthermore, D^y- 
solutions have the additional feature that they may not exist for values of a lower than crjv(Vo). 
This was explained above in detail for the case A^ = 1. Figure ^ illustrates this generic behavior 
for a D2-solution. 

As for A = 1, for this family there are generically two main types of potential well, depending 
whether ui^q is larger (shallow potential) or smaller (deep potential) than rimin (left and right parts 
of Figure |8|, respectively). In the case of a shallow potential the D2 solution exists for any width 
below a3(Vb); this is not the case for deep potentials. We will not comment on Figure ^ in great 
detail, but we note that even in the simple case of a shallow potential interesting bifurcations occur. 
Let us focus on this case. For simplicity we discard from the discussion solutions forming cnoidal 
waves up-stream (i.e., for x > a). Then, for small widths, the only possible D2 solutions are even. 
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Figure 6: The Di solution in a strongly attractive square well of width just above ct*(Po)- Left and right 
parts as in Fig. ||. The density outside the well is composed of two half solitons. This plot corresponds to 
Vq^'^ — 0.5 and a ~ 2.4^ in a low density beam at velocity Uoo/coo = 0.7. The density evolves from rtoo (far 
down-stream) to n(G) — rimin in the potential W{n) with classical energy i?^ = W{noo)- It then evolves in 
Wq from n(0) to nmin,o and back. Finally the up-stream density oscillates in the classical potential W{n) 
from n{a) = rimin back to Uoo- 

For a certain width the minima of density occur exactly at x = and x = cr (it is easy to see that 
this width coincides with al{Vo) defined in Eq. (^T[)). From there on, the previous even solutions 
still exist, but new solutions appear. They correspond to a portion of a soliton with its minimum 
before the well and one portion of cnoidal oscillation inside the well (see Figure ^, left part). This 
solution is degenerate in the sense that there exists a symmetrical equivalent solution (where the 
minimum of the soliton occurs beyond the well). It disappears when a = o"2(Vo). For a just below 
this value, one has exactly one soliton out of the well and one cnoidal oscillation inside the well; 
hence the trough of the solitonic part of the solution is sent to infinity (a feature that is not clearly 
seen on Fig. |8| due to numerical difficulties) . For larger values of a one has to oscillate more than 
once in the region of the well, but these are D3 and not D2 solutions. On the other hand, the even 
solutions still exist until a = ^^{Vq). The situation is slightly more complicated for deep potentials, 
but the basic ingredients are the same as for shallow potentials and we present the different allowed 
density profiles in Fig. ^ (right part) without detailed discussion. 

5.2 high beam velocity: v^/coo > 1 

We now consider beam velocities Voo larger (in absolute value) than Cqo = c(noo)- In the language 
of Fig. H, n^o is in this case of type ni (the minimum of the potential) and E~i = W{ni). The 
only possible flat solution far in the rear of the obstacle is a flat and constant density, namely 
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Figure 7: Synoptic diagram of the evolution of the morphology of the Di solution for a given Vq, as cr 
increases. The insets display the density profiles n{x)/ncx)- Hatched regions indicate values of a for which 
the solution does not exist. For shallow potentials (for which rii.o > rimin, upper part of the plot) the Di 
solution exists for any width cr between and 172(^0). For deep potentials (ni.o < n„iim lower part) it 
exists only if ct € [al{Vo), (J2{Vo)]- In both cases (shallow or deep potential) at a — <72{Vo) the Di solution 
disappears and becomes a B2 solution. 

n(x < 0) = noo- The matching condition ( p6| ) yields E'^ = E~i + Vb(noo — n{a)). One should verify 
that > W{ni = noo) = (see Fig. ||) and this imposes n{a) < n^. This inequality is trivially 
satisfied, because the solution inside the well is a cnoidal Wcive with. TiyyiQ^x — ^oo 

(see Fig. g). If 

we denote by 712 (as in Sec. ^) the second zero of dW/dn (the first one being rii = rioo) one has 
also to verify that E"^ < W{n2)- This imposes 

n(a) > n^^inf = - ^ [W{n2) - VF(ni)] . (32) 

Once this condition is fulfilled, the up-stream solution is a cnoidal wave. Because these solutions 
have a constant down-stream density, we denote them C-solutions. 

Two different cases are to be considered. The first and simpler one corresponds to nmin,o > 
na,inf- It occurs for high beam velocity (when W{n2) — W{ni) is large) and for shallow potentials 
(when Vq is small). In that case, the density of the cnoidal wave inside the well oscillates between 
iT-minfi Sind n^axfli ^nd the matching at x = cr is always possible. The incoming wave (for x > o") is 
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Figure 8: Synoptic diagram of the evolution of the morphology of the D2 solution for given Vq, as a increases 
(see caption of Fig. 0). For a larger than (J3{Vq) the D2 solution disappears (it becomes a B3 solution). The 
left part corresponds to a shallow well and the right part to a deep well. For legibility, the cases denoted as 
"typical" have been taken to consist ahead of the obstacle (for a; > cr) of a portion of a soliton (whereas the 
most generic solutions are cnoidal waves). 



a cnoidal wave corresponding to oscillations of the "classical particle" in the potential W{n). This 
behavior is illustrated in Fig. |^(a). 

The other case corresponds to n^^^inf > nmin,o- Then all the widths are not acceptable, namely, 
a should be such that n{a) > ria-^mf- If one denotes by La the length the density takes to go from 
the value n^^inf to Umaxfi and by Lq the period of the cnoidal wave in the well, one has 



^'^"72 



max,0 



dA 



and 



dA 



(33) 



The allowed values of a are in the intervals [0, Lq] U [Lq — Lq, Lq + La] U [2Lo — La, 2Lo + Lq] U .... 

The region of validity of the first case nmm,o > f^a,inf can be evaluated analytically in the low 
density regime. In this regime one obtains Wiji'}) — W(n\ = rioo) = Uoc F{voo/coo) 7(4'^^) where 



Fiz) 



1 



+ 1 



3z^ 



8 

1 + ^ 



(34) 



This yields no-,m//''^oo = 1 — F {z) / [^Vq S^"^) (we set z = Woo/coo)- One also obtains nmin,o/n'c 
G{z,Vo^'^) where 



Giz,Voe 



z2 + l 



(35) 
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Figure 9: A C-solution in an attractive square well. Left and right parts as in Fig. |^. Both density plots 
on the left part concern a low density beam, with Uoo/coo — 1-6 and a square well of width a/£, = 4.68174. 
Plot (a) is the generic case and corresponds to Vb = 0.1. For these values of Vq and Vao the condition (32) 
is fulfilled for any a. Plot (b) corresponds to a well of depth Vb^^ = 0.5. The parameters have been chosen 
in this case such that E^^ = W{n2) and the C-solution is just about to disappear. The right plot illustrates 
the behavior of the typical solution (corresponding to part (a)) in the diagram {n,W{n)). The density is 
constant and equal to Uoo for a; < in W{n). From x = to x = a the "classical particle" evolves in Wo{ti) 
(with an energy E^^). Finally the up-stream density oscillates in the classical potential W{n) as a cnoidal 
wave. 



The region of validity of the condition n^^inf < nmin 
can be displayed in a diagram {z = Voo/coo, Vq ^^). 



(where all the potential widths are acceptable) 
It corresponds to the region 



1 



F(z) 



i.e., to the domain below the solid line in Fig. IT 



<Giz,Vof 



(36) 



6 The repulsive square well 

In this section we consider a simple repulsive potential, namely, we take V|| (x) to be zero everywhere 
except in a finite region < x < a where it takes the constant value Vq (Vq > 0). This type of 
potential cannot correspond to a bend in the guide, but it can be realized with a (blue detuned) 
far off-resonant laser field. 

Eq. ( [25| ) still holds after changing the sign of Vq. Eq. ( [2^ ) holds also, but we rewrite it here 
with the appropriate sign for future reference. 



E- - Vo n(0) = = - Vo n{a) 



(37) 
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Figure 10: Representation (for low density beams) of the part of the plane {voo, Vq) where the C solution is 
allowed for any potential width. The allowed region is located below the solid curve. We also display on the 



same diagram the region where C solutions are allowed for any width of a repulsive potential (see Sec. 6.1). 
It corresponds to the domain below the dashed curve for low density beams and below the dot-dashed one 
for high density beams. 



and similarly one has here (instead of (|27D) 

Wo{n) = W{n)-Von 



(38) 



The potential Wo(n) has, for Vq relatively small, the same behavior as W{n), namely one minimum 
at ni^o and one maximum at n2fl, with ni < riifi < uq < n2,o < n2 (the notations are defined in 
Sec. see Fig. |l|). On the other hand, for large Vq, Wq is a monotonically decreasing function of n. 
The transition between the two regimes occurs when fj, — Vo = e(no) + /{2nQ). Hence, the terms 
"low" or "large" Vq we just defined (we also speak below of "weak" and "strong" potentials) are 
not intrinsic properties of the potential, but depend on the chemical potential ^ and on the flux ^ 



of the incoming beam (this remark is made quantitative in Sec. 6.1 below). 

Amongst all the possible solutions, we choose again to select those corresponding to a flat 
density at left infinity, (n(x) Uoo as x — > — oo), with a negative velocity {v{x) — > Voo < 0). Under 
these conditions, the transmission modes through repulsive potentials appear to be simpler than 
for the attractive ones discussed in Sec. ^, as we are now going to show. 

6.1 high beam velocity: foo/coo > 1 

The stationary solutions of a supersonic beam encountering a repulsive square well have important 
similarities with the C-solutions of Sec. ^ and will be given the same name. For weak potentials 
(Vq < fJ- — e(f^o) + ^^/(2?^o)) solutions exist whatever the value of a, and their shape is very 
similar to the transmission mode illustrated in Fig. 0(a). The main difference is that here, one has 
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n{x) > n{0) = rioo in the region of the well (whereas the reverse inequality holds for C-solutions in 
an attractive well). 

On the other hand, for strong potentials the solution in the region of the well is not (as for weak 
potentials) a cnoidal wave. In this region the density increases monotonically from x = up to 
X = a; then "classical" motion occurs in W{n) and the up-stream density profile further oscillates 
as a cnoidal wave. This means that there exists a maximum value of a above which no solution 
can be found in a strong potential. When a reaches this maximum value, the solution for x > a 
is a portion of a soliton and the whole solution has a step-like shape going from = ni (far in 
the rear of the obstacle) up to n2 at x — > +00. The behavior in this limiting case is illustrated in 
Fig. |l^. There is here a difference with the case of a strong attractive square well: when a has 
exceeded this maximum value, no other stationary solution appears for larger widths. The reason 
is that in an attractive square well, the solution in the region of the well is periodic, whereas it 
increases here without limit for large wells. 

Due to the similarities with the case of an attractive potential, it is interesting here also to 
determine more precisely for which beams a stationary supersonic solution can be found (for weak 
potentials) for all values of a. An analytical treatment is here possible in both the low and high 
density regimes. The two cases can be treated on the same footing by introducing an index v with 
I' = 1 in the dilute regime and = 1/2 for high densities. One has 
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where z = foo/coo- One then obtains UQ/n^ = z^+^ (no is defined in Sec. Fig. ||). The 
condition for a repulsive potential to be considered as weak is thus — Vq > e{no) + <^^/(2no) = 
eo + {no/riooY (1 + f )/(2^2)_ ^^^^^ ^ ^ ^2^/2 + e(noo), this reads 

9 1 V z"^ 1/ JJL. , , 



2 4 2 V 2 



This region of weak potential corresponds in Fig. IC to the domain below the dashed curve in the 



dilute regime v = and below the dot-dashed curve in the high density regime v = 1/2. As stated 
in the beginning of this Section, the terms weak or strong potential do not characterize an intrinsic 
property of the well. At Voo = Cqo for instance, all the potentials are "strong" . 

In view of future experimental studies of the system, it is also interesting to determine (in the 
regime of "strong potentials") the critical value of a after which no C-solution can be found. We 



denote this value by L^. Following a reasoning similar to the one of Sec. |5.2| one obtains 



Lr = —^l = where n„,sup = ^^oo + 77 [W{n2) - W{ni)\ . (41) 

V2 J>™„.o /rO _ lyg(^) Vo 



(42) 



In the low density regime L^- can be expressed in terms of an elliptic integral: 

Lr f^(^y(^yoe) dx i [f{z)V'^ 



i Jo V2^[x2 + x(4yoe2 + l-^2)+4Vb^2]i/2 V 8 
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Figure 11: A C-solution in a repulsive square well. The parameters have been chosen in this case such that 
E'^i = W{n2) and the C-solution is just about to disappear. Left and right parts as in Fig. |. The left plot 
concerns a low density beam, with Uoo/coo — 1.6 and a square well of width a/£_ — 3 and depth Vg — 0.1. 
The right plot illustrates the behavior of this solution in the diagram {n, W{n)). The down-stream density 
is constant and equal to riao- From x = to a; = tr the "classical particle" evolves in Wo(n) (with an energy 
E^i), then it evolves back in W{n), just on the separatrix. 



Note of course, that Eq. ( |4^ ) is meaningful only when the condition ( |40| ) is violated. La in Eq. (33) 
can also be defined with a similar expression in the low density regime. Both expressions are well 
approximated by the left part of (42), meaning that one has 8(L^/aVbO^ ^ F'i^)- We will see 



below that this corresponds to approximating the potential by a ^-function (cf. Sec. 0, Eq. (|5 
6.2 low beam velocity: foo/coo < 1 

It is easy to check that in this case, only Di solutions can be observed (there are no other D- 
solutions, and B-solutions are forbidden). The down-stream solution starts at left infinity from a 
density noo which is of type 77-2 in the terminology of Sec. ^ n{x) decreases from this value, and 
in the diagram {n,W{n)) the fictitious classical particle evolves in the potential Wo{n) during a 
"time" a, and then evolves in W{n) again. Since one should verify that i?^ < W{n2) = E~^, from 
( |37| ) this imposes n{a) < n(0). If the inequality is strict, the solution for x > cr is a cnoidal wave. 
If n{a) = n(0) the final solution is a portion of a soliton. 

Let us consider a strong potential first. If the well is narrow (a — > 0), there are two possible 
solutions, depending whether, in the diagram (n,W{n)), one "jumps" rapidly or not from W to 



Wq. This is illustrated on Fig. |12[ When a increases, these two solutions merge and disappear. 
This behavior was already observed by Hakim in the case of a model repulsive (^-potential pO[ |. 
Hakim showed (for the one dimensional Gross-Pitaevskii equation, i.e., in the dilute regime) that 
the solution which "jumps late" to Wq (right part of Fig. [T^ ) is unstable, and argued convincingly 
that the other is stable. The same situation is expected to occur here. 

In the case of a weak potential there are also two types of solution, but they do not disappear 
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Figure 12: Schematic representation in the (n, W{n)) diagram of the behavior of Di-solutions in a strong 
repulsive square weU (woo/coo < !)• The solid curve represents W{n) and the dashed one Wo{n). For this 
type of potential there are two possible solutions for a given value of the width a (as discussed in the text). 
This is illustrated by the figure, namely the "time" the fictitious particle spends in the potential Wq is the 
same in the left and right plots. For simplicity, the examples we give here are drawn in the particular case 
n(0) = n{a). 



when a increases. The point is that one can have here arbitrarily close to the separatrix energy 
Wo(n2,o)) and the period of motion in the potential Wq can thus be made as large as desired. From 
(|^, one has = Wo{n2,o) if re(0) = with 

n*o = n2,o + ^ [Win2) - W{n2,o)] . (43) 

The expression "jumping soon" (or "late") from W to Wq that we used in the discussion of strong 
potentials refers, for a weak potential, to the case where n(0) > tiq (or n(0) < Uq). If n(0) > Uq, 
the density remains larger than 77-2,0 and reaches this value at its minimum in the limit a ^ 00. If 
n(0) < Uq, > Wo{n2fi) and the fictitious classical particle evolves above the separatrix in the 
potential Wq. As in the case of a strong potential, this solution is expected to be unstable. The 
other one is certainly stable since one can show that it is identical to the result of perturbation 
theory in the limit of a very weak potential. 



7 Simple solutions in the presence of an obstacle 

The aim of this section is to study, by means of perturbation theory, some simple solutions of 
Eq. ([TtI) valid for a generic potential Vj|(a;). We will argue that near the speed of sound this 
approach fails, and that in this regime any potential can be approximated by a 5-peak. We will 
then study the scattering modes of the condensate in the presence of this potential. It allows for a 
qualitative and simple understanding of the solutions obtained for more realistic potentials in the 



previous sections. Some of the results presented here have already been obtained by Hakim |3C], 
who considered repulsive potentials only, in a slightly less general setting. 
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We again restrict the analysis to those transmission modes tending to a flat density at x — > — oo, 
n{x —oo) — > Tioo, with a negative velocity Voo- These are of the form A(x) = A^o + 5A{x) (with 
^oo = '^oo)- Denoting by Cqo the sound velocity at density rioo (Eq. (pO|)), a perturbative treatment 
of Eq. (0) yields 

6 A" + 4{vl-cl)6A = 2A^V^i{x) . (44) 
The solutions of (|43) that tend to zero when x ^ — oo are of the form 



I^^y\\{y)ew{-k\x -y\}dy when foo/coo < 1 , 
5A{x) = { (45) 
/-oo ^1 (y) sin{fc(x - y)} dy when foo/coo > 1 , 

where k = 2 - c^l"*"/^. 

Denoting by a the typical range of the potential V|j , in the limit ka ^ 1 the Green function of 
(H4|) is almost a delta peak and Eqs. (Hsl) take the simple form 



5A{x) = sgn (t;oo/coo - 1) ^ ^| (x) . (46) 

This result may seem unnatural at large velocities. Indeed, for a repulsive potential for instance, 
the density increases in the region of the potential. This kind of behavior was already found in 
Ref. |31], where a very special potential was designed for which this phenomenon occurs at any 
Voo/ Coo > 1- From (|4^ ) we see that similar motional dressed states exist for any potential in the 
limit of very large velocities (when ka ^ 1). 

In the case foo/coo > 1, an asymptotic evaluation of (|45|) yields far ahead of the obstacle (in 
the limit /cx ^ 1) an amplitude of the form 

M(x) = ^Im{e^^^Vl|(fc)} + ^^|(x) + ^ (p) , (47) 

where V\\{k) = dx eyip{—ikx)V\^{x) is the Fourier transform of Vj|(x). This shows that, for 
Voo/ Coo > 1, Eq. ( ^5D corresponds to a C-solution. The fact that the up-stream solution oscillates 
with a wave-vector k corresponds to the stationarity condition Vp{k) = —Voo (see the end of Sec. ^). 
The wake is characterized by a wavelength 2-11 /k that depends on the velocity of the beam (it 
decreases when \voo\ increases). This is due to the particular form of the dispersion relation (|2^), 
and does not occur above the Landau critical velocity in liquid helium for instance, where the 
location of the roton minimum flxes the wavelength of the wake ||3^ . 

When the beam velocity is lower than the speed of sound, for attractive (resp. repulsive) 
potentials Eq. ( ^5|) describes a Bi (resp. Di) solution. In the attractive case for instance, the bump 
density measured with respect to the constant density Hoo contains a number AA^ of atoms given 
by AA'' = 6n{x) dx — (4noo/fe^) J-^ dxV\\{x). This formula diverges when foo/coo —>■ 1 
since A; — > 0. In that limit, however, the perturbative treatment is not justified. Indeed, Eq. (45) 
gives a sensible result only if \6A\ <C ^oo- Denoting by Vq the typical value of the potential V||, this 
reads 

Vq a/k when ka <^ 1 , 



Vo/k'^ < 1 when A;cr > 1 . 



(48) 
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These conditions are satisfied only if the beam velocity v^o and the sound velocity Cqo are not too 
close. 

From Eqs. (|4^) we see that for Voo/coo < 1 the typical length scale of variation of 6 A is 
proportional to k^^. If \voc\ approaches |coo| from below, this length scale diverges and the spatial 
extension of 5A increases indefinitely. In this case it is legitimate to approximate V|j by a 6- 
function. Hence, we do not pursue the perturbative treatment any longer and turn now to the 
study of solutions in a 6-peak potential. 

Consider a potential of the form Vj| (x) = A 6{x) with A positive or negative. A realistic V|| can 
be approximated by such a potential if its typical length scale a is much lower than the healing 
length ^. In this case the approximation is valid for any beam velocity. As discussed above, any 
potential V|| can be approximated by a 5- function when Voo approaches Cqo, i.e., in the limit ka <^ 1. 

The allowed transmission modes of Eq. ( p!7| ) when V|| is a (5-potential are very simple. They 
are obtained by joining together the solutions of two straight guides (see Sec. Q), one down-stream 
(x < 0) and the other up-stream (x > 0), with the matching condition 

A'{0+) - A'{0-) = 2XA{0) . (49) 

The integration constant Ed of Eq. (|l^) changes discontinuously from the value E^^ for x < to 
the value for x > (S^ = i[A'(0=^)]2 + W{n{0))). Using (||) these two values are related 
through E+ = E- +XA{0) [^'(0") + ^'(0+)]. 

2-5 I ^ 7 1 




Figure 13: Domains of existence of the different solutions occurring for a dilute beam in a potential 
V||(x) — XS{x). Region I occurs only for Voa/cao > 1 and corresponds to C-solutions. Region II corresponds 
to Bi-solutions and occupies the domain under the dashed line (it occurs only for attractive potentials and 
for Voa/coo < !)• Region III corresponds to D-solutions: Di-solutions for A > and D2 solutions when A < 0. 



Let us first consider the case foo/coo > 1- We know from Sees. [5.2| and 6.1 that one should 
observe C-solutions. One has in this case a constant down-stream density = A"^ and -E^ = 
W{noo)- The matching condition imposes E+ = E- +2X'^noo and ^'(0+) = 2Ayloo- In order to 
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have a finite solution at x = +00 one should moreover satisfy W{noo) < i^"^ < W{n2). This fixes 
an upper bound for the intensity of the perturbation, given by 2A^noo < W{n2) — W{noo)- This 
is the criterion for observing C-solutions. In the low-density limit this relation takes the analytical 
form 



;X^f<F{^] , (50) 

Coo 



where the function F is defined in (|34|). It corresponds to the region I of Fig. |13. 

Let us now consider the case foo/coo < 1- The first and simpler solution is the Bi-solution 
which occurs for attractive potentials. The corresponding matching condition ( [49D reads in this 
case n'(0)/n(0) = —2 A. It can be fulfilled whatever the value of A (A < 0) because the function 
(n'/n)^ = 8n~^[Eci — W{n)] can be made arbitrarily large (it grows as 8e{n)/n for large n). The 
domain of existence of Bi-solutions corresponds to region II in Fig. 13. 

D-solutions can also be observed if t'oo/coo < 1- For A > these are the Di solutions whose 
behavior in the {n,W{n)) diagram is illustrated in Fig. |l^ (in the case of a depth of finite width 
though). For A < these are D2-solutions such as those shown in Fig. ^ for a ^ 0. The limiting 
case (beyond which these solutions disappear) is obtained for two symmetric portions of solitary 
waves. The matching condition (^) then fixes the maximum value of |A|, which corresponds in 
both cases to 

\X\^ = k(^) where K{z) = ^ \-8 - 20 + 1 + {I + 8 z^/^V^^ (51) 

V Coo / 4: Z ^ J 

The corresponding domain is shown in Fig. |l^ (region III). 



8 Experimental considerations 

In view of future experimental observation of the different flows past an obstacle, we now evaluate 
the orders of magnitude of the relevant dimensionless parameters identified above. For concreteness 
we consider a beam such as the one in preparation at the ENS |Q : ^^Rb atoms are guided along 
the X direction, with a harmonic transverse confinement {uj± = 2 vr x 500 Hz and a± = 0.5 ^m). The 
beam has a velocity Voo of order of 0.5 m/s and a flux $ varying from 10^ to 10^ atom.s"^. Hence the 
quantity rioo flsc varies from 10~^ to 1. For a rough estimate of the order of magnitude of the relevant 
parameters, we will consider that this corresponds to the low density limit rioo o-sc ^ 1- Then, the 
healing length (defined in (^)) and the speed of sound vary from ^ = 25 ;um and Coo = 20^m.s^^ 
(for noo isc = 10~^) to ^ = 0.25 /im and Coo = 2 mm.s""*^ (tIoo a^c = 1)- Note that the beam velocity 
is much larger than the speed of sound (the quantity Woo/coo is of order 10^ at least). 

If the obstacle is a bend of constant radius of curvature Rc and opening angle 0, Vq = —1/(8 i?^) 
and a = @ Rc- A reasonable order of magnitude is Rc = 5 a±, leading to Vq^^ = 2S0 (^/'^-'-)^ which 
varies roughly from 10~^ to 10. For O = Tr/2 one has a = 8a± and the obstacle can safely be 
treated perturbatively because in this case ka ~ 10^ ^ 1 and Eq. (^) is satisfied since 

Vo/k'^ « 10-^ < 1 . 

For this configuration one thus expects profiles in agreement with Eq. (|45| ) (case Uoo/coo > 1), i-e. 
C-solutions. Also, since we are in the regime ka ^ 1 Eq. ( ^6[ ) holds, meaning that the wake ahead 
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of the obstacle is very weak, and that there is a decrease in the density in the region of the bend. 
However, this decrease is extremely small: it corresponds to a number of atoms smaller than unity, 
and under these conditions nothing noticeable is expected to occur in the bend. 

The situation changes drastically if the obstacle is due to a transverse laser beam because the 
potential can be made much stronger. It can moreover be attractive or repulsive depending on 
the laser's frequency, and also the velocity of the obstacle relative to the beam can be modified by 
using an acousto-optic deflector. 

We consider a laser with power P = 70 mW and a wave length varying from 780 to 790 
nm (the atomic transition corresponds to a wavelength = 780.2 nm and has a natural width 
r = 12 7r MHz). The laser beam has a typical waist a of the order of 50 fim. Then one obtains for 
the transverse potential 

V\\{x) = VoeM-^} where Vq = f, , (52) 

S = iOji — u!^ being the detuning and Is = 16 W.m"^ is the saturation intensity. This yields 

In this case one is in the regime ka ^ I, with Vo/fc^ ~ 10"^ or 10~^. Hence one is on the edge of 
applicability of the perturbative approach. By changing the wavelength of the laser or its velocity 
with respect to the beam, one may enter into the non-perturbative regime. One can still find C- 
solutions for this type of potential if the parameters of the system remain in the appropriate region 
of figure C-solutions exist for any value of cr if Vq < ]^ (^oo/coo)^ in the case of an attractive 
potential, and if Vq < 5 (t'oo/coo)^ in the case of repulsive potential (since foo/coo 1 we consider 
here the asymptotic versions of Eqs. ( |36| ) and (|40|)). By modifying the value of the healing length 
and of the beam velocity these conditions can easily be satisfied or violated (especially in the case 
of repulsion). In the region where these conditions are violated one can experimentally study the 
transition from a stationary flow (of type C) to a time dependent one. For instance, it would be of 
great interest to study the modification of the drag at the transition. 

Just at the boundary between the two regions, for a repulsive potential the density profile 
has the behavior shown in Fig. |ll| (for an attractive potential the density profile corresponds to 
Fig. |^(b)). The density has a step-like shape, and the beam ahead of the obstacle has a velocity 
lower than the speed of sound. One has just at the transition a = (or La in the case of an 
attractive potential) with 8{VoLr^)'^ ^ F{voo/coo) ~ (woo/coo)^/4 (see the end of Sec. 3.1). From 



the estimate ( |53|) this occurs for i/cr ^ 10^^ or 10"'^. For a laser, one can tune the waist a by a 
factor of order 5 say, and more important, by changing the density one can modify the value of 
and indeed reach the appropriate regime. This would have a very important effect on the beam, 
since the velocity ahead of the obstacle would be lower than the speed of sound whereas it is of the 
order of 1 m.s~^ down-stream. Accordingly, the density along the beam would go from = n2 
to ni (see Fig. ^l]). It is not difficult to see that, in the dilute regime, in the limit Uoo/coo ^ 1 one 
has 77,2 — n,i(t;oo/coo)^/4. Thus the down-stream beam density is divided by a factor of order 10^ 
with respect to the up-stream one and, by conservation of the flux, the velocity is multiplied by 
the same factor. The beam velocity ahead of the obstacle is then of the order of a few micron per 
second! 
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9 Conclusions 



In this paper we have studied the different stationary profiles of a Bose-Einstein condensed beam 
propagating through a guide with an obstacle. The beam far down-stream is characterized by its 
velocity Voo and density Uoo (or equivalently, its healing length ^). The obstacle is represented by 
a one dimensional square well of depth itVo and width a. Let us denote by A the region along 
the guide where the potential is different from zero. Our study allows to identify the relevant 
dimensionless parameters governing the flow across the obstacle. These are n^agc, Vao/coo, a/^ 
and Vq . By varying these parameters one obtains a wide range of density profiles (we identify 
three main stationary families denoted as B, C and D-solutions). We have numerically checked 
that similar results are obtained for potentials other than square wells (such as V||(3;) = 
zbVo exp{— x^/o"^} for instance) and we thus believe that our analysis of the flow is quite general. 

The richest case occurs when the external potential is attractive. In the subsonic regime (beam 
velocity lower than the corresponding speed of sound) the simplest solution is a symmetric density 
having a peak in A (B-family). These solutions may have density oscillations in A, and are always 
symmetric. Another type of transmission mode is a soliton-like depressed solution pinned to the 
obstacle (D-family), which may also have density oscillations in A and, unlike the B-family, a wake 
up-stream. Finally, the supersonic transmission modes (C-solutions) possess, in the simplest case, 
a density trough in A and are constant outside. They may also have density oscillations in A and 
an up-stream wake. 

For a repulsive potential, in the subsonic regime the transmission modes are of type D with no 
density oscillations in A. In the supersonic case, the modes are of type C but with a density peak 
instead of a trough. Step-like solutions of increasing density across the obstacle (with or without 
wake) also exist. Specifically, we have identified an interesting (and experimentally reachable) 
regime where the beam is almost stopped by a repulsive obstacle and gains ahead of it several 
orders of magnitude in density (see Sec. ^). 

An important aspect of the problem that remains open are considerations related to the stability 
of the solutions. Some related work is in progress l33]. We just note here that the limiting C-solution 



shown in Figs, ^b and |11| can be turned into a D-solution by exchanging the down-stream and the 
up-stream behavior of the fiow. However, these are part of a continuous family of fiow patterns and 
are probably instable. The C-solutions on the other hand, as selected by the radiation condition 
(see the end of Sec. are the only acceptable ones for foo/coo > 1- Another aspect concerns Bose 
gases with attractive intra-atomic interactions, which may also be treated with our formalism. 
In this case the potential W{n) introduced in Sec. ^ has a single- well shape, and the number of 
different transmission modes whose density tends to a constant at the input of the guide is greatly 
reduced with respect to the case of repulsive interactions considered here. 

We conclude by noting that a branch of BEC that is now expanding is the non-linear counter- 
part of transport experiments of mesoscopic physics in condensed matter. In the latter case the 
coherent transport of two-dimensional electron gases through various geometries has been consid- 
ered in great detail. Future developments of transport experiments of Bose condensates should 
extend those to non-linear regimes. 
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Appendix A 

In this Appendix we recall the conditions for existence of a bound state in a bend of an ordinary 
wave guide (i.e. without non-linear terms in the Schrodinger equation) and put the adiabatic limit 
used in the text on firmer mathematical grounds. 

Let us consider a curve C, of parametric equation fc{x), x being a curvilinear abscissa along 
C. The Frenet frame {t,n,b), the curvature k{x) and the torsion t(x) are defined by t = dfc/dx, 
dt/dx = Kfi, dfi/dx = —Kt + r 6 and db/dx = — rn. 

We first introduce a curvilinear coordinate system. The position of a point of space is specified 
by coordinates (x, y, z) through 



f{x,y, z) = fc{x) + y N + z B , 



(Al) 



where N{x) = cosO n+sinO b, B{x) = — sin0 n+cos 9 b and 9{x) is defined through dO /dx = —t{x). 

We then select a potential of the form V^_{y,z). Note that the choice of vectors N and B for 
defining the transverse coordinates y and z is not irrelevant, i.e., the manner in which V±_ winds 
round C does matter: it has to be the same as the way [N ^ B) winds around t. Indeed one can 
see that some torsion may create a repulsive potential along x which could cancel the localizing 
effect of the bend. This is avoided with the type of coordinate dependence we have chosen. A 
simple way of seeing this is by noticing that V has coordinates {h~^dx,dy,dz) in the {t,N,B) 
frame (with h{x, y,z) = 1 — k (y cos 9 — z cos 9)) and the force —'W±{y, z) thus has no tangential 
component. It would have been more natural to define the (y, z) coordinates as in ( [Al| ) but using 
the {n,b) vectors instead of {N,B). In this case however the force —W±{y,z) would generically 
have a tangential component which could spoil the localizing properties of the bend. Note that the 
latter discussion is of course irrelevant for a potential V± with circular symmetry (as any simple 
experimental wave guide is expected to have). 

Denoting h' = dxh and h" = d^^h, the Schrodinger equation for <I> = /i^/^ ^ reads 



1^2 ~ yy ^ ■ 



5{h 



/\2 



8/i2 



+ ^ + V^iy,z) 



<^ = l_L<^> . (A2) 



Note that in (A2) as in the rest of the paper, we take h = m = 1. The choice of coordinates ( |A1| ) 
gives a volume element d^r = hdxdydz and ^ is thus normalized as J dxdydz]^]"^ = 1. 

The adiabatic limit is defined by/i^l,K^/i',K^ %//;/'. In this limit Eq. ( [A2| ) decouples 
into a longitudinal and a transverse equation. One obtains for the longitudinal equation a potential 
V||(a;) = —k'^{x)/8, attractive in the region of the bend. Since any attractive potential in one 



dimension has a bound state |22|, in this limit there exists a quantum state localized in the bend. 



The theorem of Goldstone and Jaffe |19| establishes the existence of a bound state for much 
more general wave-guides, with arbitrary curvature. It was originally demonstrated for sharp-walls 



30 



wave guides, but the proof can be straightforwardly extended to the case of a smooth confining 
potential. Hence we do not reproduce its derivation here. We just note that for some points of 
space, the coordinate system ( [Al| ) can be ambiguous. This imposes for the sharp wall problem 
the requirement that the transverse size of the guide does not exceed the radius of curvature of C. 
For the smooth potential problem, the same restriction requires that the typical range of V± (or 
equivalently the spatial extension of (f)±) is lower than the radius of curvature of C. Obviously this 
condition is much weaker than the condition of adiabaticity. Note also that the Goldstone-Jaffe 
theorem has the same limitation as above, namely some torsion may destroy the localizing effect 
of the bend and the way V± winds around the curve C is not irrelevant. A potential V±{y, z) where 
the coordinates y and z are defined above (Eqs. (|A1|)) ensures the applicability of the theorem. 
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